The HLA region in ANCA-associated vasculitis: characterisation of genetic associations in a Scandinavian patient population.

OBJECTIVE: The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitides (AAV) are inflammatory disorders with ANCA autoantibodies recognising either proteinase 3 (PR3-AAV) or myeloperoxidase (MPO-AAV). PR3-AAV and MPO-AAV have been associated with distinct loci in the human leucocyte antigen (HLA) region. While the association between MPO-AAV and HLA has been well characterised in East Asian populations where MPO-AAV is more common, studies in populations of European descent are limited. The aim of this study was to thoroughly characterise associations to the HLA region in Scandinavian patients with PR3-AAV as well as MPO-AAV. METHODS: Genotypes of single-nucleotide polymorphisms (SNPs) located in the HLA region were extracted from a targeted exome-sequencing dataset comprising Scandinavian AAV cases and controls. Classical HLA alleles were called using xHLA. After quality control, association analyses were performed of a joint SNP/classical HLA allele dataset for cases with PR3-AAV (n=411) and MPO-AAV (n=162) versus controls (n=1595). Disease-associated genetic variants were analysed for association with organ involvement, age at diagnosis and relapse, respectively. RESULTS: PR3-AAV was significantly associated with both HLA-DPB1*04:01 and rs1042335 at the HLA-DPB1 locus, also after stepwise conditional analysis. MPO-AAV was significantly associated with HLA-DRB1*04:04. Neither carriage of HLA-DPB1*04:01 alleles in PR3-AAV nor of HLA-DRB1*04:04 alleles in MPO-AAV were associated with organ involvement, age at diagnosis or relapse. CONCLUSIONS: The association to the HLA region was distinct in Scandinavian cases with MPO-AAV compared with cases of East Asian descent. In PR3-AAV, the two separate signals of association to the HLD-DPB1 region mediate potentially different functional effects.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Relation between HLA and copy number variation of steroid 21-hydroxylase in a Swedish cohort of patients with autoimmune Addison's disease.

OBJECTIVE: Autoantibodies against the adrenal enzyme 21-hydroxylase is a hallmark manifestation in autoimmune Addison's disease (AAD). Steroid 21-hydroxylase is encoded by CYP21A2, which is located in the human leucocyte antigen (HLA) region together with the highly similar pseudogene CYP21A1P. A high level of copy number variation is seen for the 2 genes, and therefore, we asked whether genetic variation of the CYP21 genes is associated with AAD. DESIGN: Case-control study on patients with AAD and healthy controls. METHODS: Using next-generation DNA sequencing, we estimated the copy number of CYP21A2 and CYP21A1P, together with HLA alleles, in 479 Swedish patients with AAD and autoantibodies against 21-hydroxylase and in 1393 healthy controls. RESULTS: With 95% of individuals carrying 2 functional 21-hydroxylase genes, no difference in CYP21A2 copy number was found when comparing patients and controls. In contrast, we discovered a lower copy number of the pseudogene CYP21A1P among AAD patients (P = 5 × 10-44), together with associations of additional nucleotide variants, in the CYP21 region. However, the strongest association was found for HLA-DQB1*02:01 (P = 9 × 10-63), which, in combination with the DRB1*04:04-DQB1*03:02 haplotype, imposed the greatest risk of AAD. CONCLUSIONS: We identified strong associations between copy number variants in the CYP21 region and risk of AAD, although these associations most likely are due to linkage disequilibrium with disease-associated HLA class II alleles.

Evolutionary constraint and innovation across hundreds of placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Bayesian mixed model analysis uncovered 21 risk loci for chronic kidney disease in boxer dogs.

Chronic kidney disease (CKD) affects 10% of the human population, with only a small fraction genetically defined. CKD is also common in dogs and has been diagnosed in nearly all breeds, but its genetic basis remains unclear. Here, we performed a Bayesian mixed model genome-wide association analysis for canine CKD in a boxer population of 117 canine cases and 137 controls, and identified 21 genetic regions associated with the disease. At the top markers from each CKD region, the cases carried an average of 20.2 risk alleles, significantly higher than controls (15.6 risk alleles). An ANOVA test showed that the 21 CKD regions together explained 57% of CKD phenotypic variation in the population. Based on whole genome sequencing data of 20 boxers, we identified 5,206 variants in LD with the top 50 BayesR markers. Following comparative analysis with human regulatory data, 17 putative regulatory variants were identified and tested with electrophoretic mobility shift assays. In total four variants, three intronic variants from the MAGI2 and GALNT18 genes, and one variant in an intergenic region on chr28, showed alternative binding ability for the risk and protective alleles in kidney cell lines. Many genes from the 21 CKD regions, RELN, MAGI2, FGFR2 and others, have been implicated in human kidney development or disease. The results from this study provide new information that may enlighten the etiology of CKD in both dogs and humans.

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis.

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated.

Interpretable machine learning identifies paediatric Systemic Lupus Erythematosus subtypes based on gene expression data.

Transcriptomic analyses are commonly used to identify differentially expressed genes between patients and controls, or within individuals across disease courses. These methods, whilst effective, cannot encompass the combinatorial effects of genes driving disease. We applied rule-based machine learning (RBML) models and rule networks (RN) to an existing paediatric Systemic Lupus Erythematosus (SLE) blood expression dataset, with the goal of developing gene networks to separate low and high disease activity (DA1 and DA3). The resultant model had an 81% accuracy to distinguish between DA1 and DA3, with unsupervised hierarchical clustering revealing additional subgroups indicative of the immune axis involved or state of disease flare. These subgroups correlated with clinical variables, suggesting that the gene sets identified may further the understanding of gene networks that act in concert to drive disease progression. This included roles for genes (i) induced by interferons (IFI35 and OTOF), (ii) key to SLE cell types (KLRB1 encoding CD161), or (iii) with roles in autophagy and NF-κB pathway responses (CKAP4). As demonstrated here, RBML approaches have the potential to reveal novel gene patterns from within a heterogeneous disease, facilitating patient clinical and therapeutic stratification.

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA.

OBJECTIVE: To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV). METHODS: Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay. RESULTS: PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele. CONCLUSION: We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

Contribution of Rare Genetic Variation to Disease Susceptibility in a Large Scandinavian Myositis Cohort.

OBJECTIVE: Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of complex autoimmune conditions characterized by inflammation in skeletal muscle and extramuscular compartments, and interferon (IFN) system activation. We undertook this study to examine the contribution of genetic variation to disease susceptibility and to identify novel avenues for research in IIMs. METHODS: Targeted DNA sequencing was used to mine coding and potentially regulatory single nucleotide variants from ~1,900 immune-related genes in a Scandinavian case-control cohort of 454 IIM patients and 1,024 healthy controls. Gene-based aggregate testing, together with rare variant- and gene-level enrichment analyses, was implemented to explore genotype-phenotype relations. RESULTS: Gene-based aggregate tests of all variants, including rare variants, identified IFI35 as a potential genetic risk locus for IIMs, suggesting a genetic signature of type I IFN pathway activation. Functional annotation of the IFI35 locus highlighted a regulatory network linked to the skeletal muscle-specific gene PTGES3L, as a potential candidate for IIM pathogenesis. Aggregate genetic associations with AGER and PSMB8 in the major histocompatibility complex locus were detected in the antisynthetase syndrome subgroup, which also showed a less marked genetic signature of the type I IFN pathway. Enrichment analyses indicated a burden of synonymous and noncoding rare variants in IIM patients, suggesting increased disease predisposition associated with these classes of rare variants. CONCLUSION: Our study suggests the contribution of rare genetic variation to disease susceptibility in IIM and specific patient subgroups, and pinpoints genetic associations consistent with previous findings by gene expression profiling. These features highlight genetic profiles that are potentially relevant to disease pathogenesis.

Association of Protective HLA-A With HLA-B∗27 Positive Ankylosing Spondylitis.

OBJECTIVES: To further elucidate the role of the MHC in ankylosing spondylitis by typing 17 genes, searching for HLA-B∗27 independent associations and assessing the impact of sex on this male biased disease. METHODS: High-confidence two-field resolution genotyping was performed on 310 cases and 2196 controls using an n-1 concordance method. Protein-coding variants were called from next-generation sequencing reads using up to four software programs and the consensus result recorded. Logistic regression tests were applied to the dataset as a whole, and also in stratified sets based on sex or HLA-B∗27 status. The amino acids driving association were also examined. RESULTS: Twenty-five HLA protein-coding variants were significantly associated to disease in the population. Three novel protective associations were found in a HLA-B∗27 positive population, HLA-A∗24:02 (OR = 0.4, CI = 0.2-0.7), and HLA-A amino acids Leu95 and Gln156. We identified a key set of seven loci that were common to both sexes, and robust to change in sample size. Stratifying by sex uncovered three novel risk variants restricted to the female population (HLA-DQA1∗04.01, -DQB1∗04:02, -DRB1∗08:01; OR = 2.4-3.1). We also uncovered a set of neutral variants in the female population, which in turn conferred strong effects in the male set, highlighting how population composition can lead to the masking of true associations. CONCLUSION: Population stratification allowed for a nuanced investigation into the tightly linked MHC region, revealing novel HLA-B∗27 signals as well as replicating previous HLA-B∗27 dependent results. This dissection of signals may help to elucidate sex biased disease predisposition and clinical progression.

The ABCC4 gene is associated with pyometra in golden retriever dogs.

Pyometra is one of the most common diseases in female dogs, presenting as purulent inflammation and bacterial infection of the uterus. On average 20% of intact female dogs are affected before 10 years of age, a proportion that varies greatly between breeds (3-66%). The clear breed predisposition suggests that genetic risk factors are involved in disease development. To identify genetic risk factors associated with the disease, we performed a genome-wide association study (GWAS) in golden retrievers, a breed with increased risk of developing pyometra (risk ratio: 3.3). We applied a mixed model approach comparing 98 cases, and 96 healthy controls and identified an associated locus on chromosome 22 (p = 1.2 × 10-6, passing Bonferroni corrected significance). This locus contained five significantly associated SNPs positioned within introns of the ATP-binding cassette transporter 4 (ABCC4) gene. This gene encodes a transmembrane transporter that is important for prostaglandin transport. Next generation sequencing and genotyping of cases and controls subsequently identified four missense SNPs within the ABCC4 gene. One missense SNP at chr22:45,893,198 (p.Met787Val) showed complete linkage disequilibrium with the associated GWAS SNPs suggesting a potential role in disease development. Another locus on chromosome 18 overlapping the TESMIN gene, is also potentially implicated in the development of the disease.

A novel canine reference genome resolves genomic architecture and uncovers transcript complexity.

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

Genetic and clinical basis for two distinct subtypes of primary Sjögren's syndrome.

OBJECTIVES: Clinical presentation of primary Sjögren's syndrome (pSS) varies considerably. A shortage of evidence-based objective markers hinders efficient drug development and most clinical trials have failed to reach primary endpoints. METHODS: We performed a multicentre study to identify patient subgroups based on clinical, immunological and genetic features. Targeted DNA sequencing of 1853 autoimmune-related loci was performed. After quality control, 918 patients with pSS, 1264 controls and 107 045 single nucleotide variants remained for analysis. Replication was performed in 177 patients with pSS and 7672 controls. RESULTS: We found strong signals of association with pSS in the HLA region. Principal component analysis of clinical data distinguished two patient subgroups defined by the presence of SSA/SSB antibodies. We observed an unprecedented high risk of pSS for an association in the HLA-DQA1 locus of odds ratio 6.10 (95% CI: 4.93, 7.54, P=2.2×10-62) in the SSA/SSB-positive subgroup, while absent in the antibody negative group. Three independent signals within the MHC were observed. The two most significant variants in MHC class I and II respectively, identified patients with a higher risk of hypergammaglobulinaemia, leukopenia, anaemia, purpura, major salivary gland swelling and lymphadenopathy. Replication confirmed the association with both MHC class I and II signals confined to SSA/SSB antibody positive pSS. CONCLUSION: Two subgroups of patients with pSS with distinct clinical manifestations can be defined by the presence or absence of SSA/SSB antibodies and genetic markers in the HLA locus. These subgroups should be considered in clinical follow-up, drug development and trial outcomes, for the benefit of both subgroups.

Targeted sequencing reveals the somatic mutation landscape in a Swedish breast cancer cohort.

Breast cancer (BC) is a genetically heterogeneous disease with high prevalence in Northern Europe. However, there has been no detailed investigation into the Scandinavian somatic landscape. Here, in a homogeneous Swedish cohort, we describe the somatic events underlying BC, leveraging a targeted next-generation sequencing approach. We designed a 20.5 Mb array targeting coding and regulatory regions of genes with a known role in BC (n = 765). The selected genes were either from human BC studies (n = 294) or from within canine mammary tumor associated regions (n = 471). A set of predominantly estrogen receptor positive tumors (ER + 85%) and their normal tissue counterparts (n = 61) were sequenced to ~ 140 × and 85 × mean target coverage, respectively. MuTect2 and VarScan2 were employed to detect single nucleotide variants (SNVs) and copy number aberrations (CNAs), while MutSigCV (SNVs) and GISTIC (CNAs) algorithms estimated the significance of recurrent somatic events. The significantly mutated genes (q ≤ 0.01) were PIK3CA (28% of patients), TP53 (21%) and CDH1 (11%). However, histone modifying genes contained the largest number of variants (KMT2C and ARID1A, together 28%). Mutations in KMT2C were mutually exclusive with PI3KCA mutations (p ≤ 0. 001) and half of these affect the formation of a functional PHD domain. The tumor suppressor CDK10 was deleted in 80% of the cohort while the oncogene MDM4 was amplified. Mutational signature analyses pointed towards APOBEC deaminase activity (COSMIC signature 2) and DNA mismatch repair (COSMIC signature 6). We noticed two significantly distinct patterns related to patient age; TP53 being more mutated in the younger group (29% vs 9% of patients) and CDH23 mutations were absent from the older group. The increased somatic mutation prevalence in the histone modifying genes KMT2C and ARID1A distinguishes the Swedish cohort from previous studies. KMT2C regulates enhancer activation and assists tumor proliferation in a hormone-rich environment, possibly pointing to a role in ER + BC, especially in older cases. Finally, age of onset appears to affect the mutational landscape suggesting that a larger age-diverse population incorporating more molecular subtypes should be studied to elucidate the underlying mechanisms.

SweHLA: the high confidence HLA typing bio-resource drawn from 1000 Swedish genomes.

There is a need to accurately call human leukocyte antigen (HLA) genes from existing short-read sequencing data, however there is no single solution that matches the gold standard of Sanger sequenced lab typing. Here we aimed to combine results from available software programs, minimizing the biases of applied algorithm and HLA reference. The result is a robust HLA population resource for the published 1000 Swedish genomes, and a framework for future HLA interrogation. HLA 2nd-field alleles were called using four imputation and inference methods for the classical eight genes (class I: HLA-A, HLA-B, HLA-C; class II: HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRB1). A high confidence population set (SweHLA) was determined using an n-1 concordance rule for class I (four software) and class II (three software) alleles. Results were compared across populations and individual programs benchmarked to SweHLA. Per gene, 875 to 988 of the 1000 samples were genotyped in SweHLA; 920 samples had at least seven loci called. While a small fraction of reference alleles were common to all software (class I = 1.9% and class II = 4.1%), this did not affect the overall call rate. Gene-level concordance was high compared to European populations (>0.83%), with COX and PGF the dominant SweHLA haplotypes. We noted that 15/18 discordant alleles (delta allele frequency >2) were previously reported as disease-associated. These differences could in part explain across-study genetic replication failures, reinforcing the need to use multiple software solutions. SweHLA demonstrates a way to use existing NGS data to generate a population resource agnostic to individual HLA software biases.

Common genetic variation in the autoimmune regulator (AIRE) locus is associated with autoimmune Addison's disease in Sweden.

Autoimmune Addison's disease (AAD) is the predominating cause of primary adrenal failure. Despite its high heritability, the rarity of disease has long made candidate-gene studies the only feasible methodology for genetic studies. Here we conducted a comprehensive reinvestigation of suggested AAD risk loci and more than 1800 candidate genes with associated regulatory elements in 479 patients with AAD and 2394 controls. Our analysis enabled us to replicate many risk variants, but several other previously suggested risk variants failed confirmation. By exploring the full set of 1800 candidate genes, we further identified common variation in the autoimmune regulator (AIRE) as a novel risk locus associated to sporadic AAD in our study. Our findings not only confirm that multiple loci are associated with disease risk, but also show to what extent the multiple risk loci jointly associate to AAD. In total, risk loci discovered to date only explain about 7% of variance in liability to AAD in our study population.